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ABSTRACT 


An exverimental study was carried out to investigate 
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Nicolet storage oscilloscone. The CHS computer system was 
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The results concluded that the temverature of the 
boundary layer gas had the controlling effect over the 
cylinder bulk gas temperature in determining the actual rate 
of heat transfer through the cylinder walls. The apparent 
zero rate of heat transfer point occurred when the bulk gas 
temperature equaled the cylinder wall temperature. Because 
Siecne strong effect that the boundary layer had on the 
heat transfer in the cylinder, tne actual zero rate of neat 
transfer volnt occurred when the boundary layer gas tenper- 
ature equaled the cylinder wall temoerature. The boundary 
layer gas temperature was determined responsible for up to a 
64 degree phase difference in the times of the actual and 
apparent zero rate of heat transfer points in the compression 
and exvansion strokes in the cylinder. This vhase difference 
was evaluated to be leading in both cases. Therefore the 
actual time of zero rate of heat transfer preceded tne an- 
parent time of zero rate of neat transfer by a significant 
margin in each piston stroke, 
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CHAPTER T 
INTRODUCTION 


The reciprocating viston-cylinder system 1s a common 
and relatively simple device used in many different apvlica- 
tions in today's mechanical world. Each day some type of | 
reciprocating piston-cylinder system impacts on our lives 
whether it be through a reciprocating automotive engine or 
through some other type of reciprocating pumping device that 
may provide another personal service. An additional example 
of this type of device is a reciprocating cryogenic expander, 
whose purpose is to produce a reduction in the temperature 
of the gas passing through it. 

In all of these tyves of systems work is vroduced 
marouesn the expansion of a tluid against a piston. For an 
ideal expansion engine the exvansion stage should be isen- 


tropic and the work per unit mass should be expressed as: 
1 = - 
Ns Cy eo) (1) 


However, because this fluid expansion is not an ideal pro- 
cess, losses and inefficiencies are evident. Some examples 
of typical loss mechanisms usually associated with through 
flow expanders are: 

i. clearance losses = Because of the clearance 
volume remaining when the yviston is at tov dead center, the 


cylinder cannot fully exhaust the entire old charge. There- 





fore the fresh charge entering the piston is further limited 
in the actual volume that it can ideally occupy. The clear~ 
ance losses then define these variations from the ideal 
casee The actual factors responsible for the clearance loss 
include much more complex processes such as blow in-blow out 
Mess ecitects and the intricate influences of the mass flow 
rate. These processes are extremely involved and their 
evaluation here is not required, so then just calling atten- 
haon tO their existence will suffice. 

ee valve losses ~ Because the inlet and exhaust 
valves are resvonsible for a resistance to the fluid flow, 
there will be a pressure drop across the valves tnat will 
account for some loss. 

3e piston ring losses = There will be some loss 
associated with the leakage vast the piston rings. Also 
any friction between the vniston rings and cylinder walls 
will generate some heat that will be transferred to the 
memid and add to the loss. 
The avove examples are certainiy not the only loss mech- 
anisms evident in a piston-cylinder arrangement, but they 
are usually considered the prominent factors. 

This experimental investigation will consider the loss 
mechanisms evident in a closed cylinder. No inlet or ex- 
naust processes will be evaluated. The mass of the working 
fluid will remain constant throughout the compression expan- 


Sion cycle being examined. Therefore the causes of the 


10 





Losses evident in the cylinder will simplify, as the neat 
erpomster loss tall be the only vrominent loss mechanism. 
This will then vermit the determination of the actual effect 


that the heat transfer loss has on the system. 


tei. Calculation of Loss 


To evaluate the thermodynamic vrocesses occurring in 
the closed piston=cylinder the first law of thermodynamics 
fete se Utllized. in its most basic form the first law 


states that 
dG = dU + al. (2) 


If the pressure of the gas is spatially uniform at each 
state of the cycle and the work done by the gas is re-= 
perictved tO the displacement of the piston, then the work 


term can be defined as 
ay = { Pav. aS, 
For an ideal gas the change in internal energy of the gas is 


dU = mC, ati Paes (4) 


The first law can then be exvoressed as 


fray = dq - me, dou oe (5) 


if the process is assumed to be ideal and adiabatic, 


then the heat transfer in the system will be zero (dQ = 0) 


is 





Bema tne first law will then simplify to 


[Pav = mC. alii eece (6) 


When dT is eliminated by substitution of tne ideal gas 


equation of state, equation (5) integrates to 
S 
PV = constant. (7) 


mae Construction of the trace for this process on a vres- 
Buse VS. Volume dlagram Will yield a plot in which the con-= 
pression and expansion processes share the same line and 
therefore describe no enclosed area (Figure 1). There will 
be zero net work done on the working fluid. 

If the process is not ideal, the compression-exvansion 
meace Will not follow the py° =—OCMStancl patie  PLOCLing 
the non-ideal expansion trace on a pressure vse volume 
diagram will yield a closed plot circumscribing some en- 
closed area (Figure 2). This enclosed area revresents the 
net work done on the working fluid during the cycle and also 
represents the total losses that the system experiences 
during the cycle. 

Therefore the total loss per cycle in the piston- 
cylinder arrangement can be expressed as the integral of 


pressure times a differential volume around the entire cycle 
f Pay = loss. (8) 
This loss is specifically associated with just the bulk gas 


ike 
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in the cylinder. The total loss ver cycle can also be ex- 
oressed in a different form by considering the loss as the 


integral of the rate of heat transfer from the bulk 


093 
§ 
U) 


through the wall around the entire cycle 


pa dt = loss. (9) 
ioe boundary Layer Sfifects 


Adjacent to the inside of the cylinder wall a small 
volume of gas with finite mass will permanently exist 
Mamoush whitch the heat flux must travel in order to reach 
meme cylinder wall. As the viston moves in the cylinder, the 
G@epth Of this thin layer of gas will vary constantly due to 
Mell defined coundary layer effects. If the cyclic heat 
flow entering and exiting the boundary layer is considered, 
the importance of the ocoundary layer effect can be illus-~ 
Peeced (Figure 3). 

For an entire cycle the total loss from the bulk gas to 
the wall (pa dt) should be equal to the loss from the bulk 
gas to the boundary layer ga, dt) and should also be eauiv- 


alent to the loss from the boundary laver to the wall 


(Paz aoe 
ga dt = pa, dt = $45 dt (10) 


Because the heat capacitance of the cylinder wall is very 
great when compared to the heat capacitance of the gas, 


the wall temperature can be assumed to be constant when 
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considering a single cycle. 

If the small volume of gas comvorising the boundary 
layer had no heat capacitance then at any instant during the 
cycle the rate of heat transfer from the bulk gas to the 
boundary layer would equal the rate of heat transfer from — 


the boundary layer to the cylinder wall. 


are agence tt zero heat 2 11) 


capacitance in boundary laye 

However, the presence of this boundary layer would suggest 
Maat it would hold some finite value of heat capacitance 
ovetween the temperature extremes of the bulx gas and the 
cylinder wall. Assuming that the doundary layer does suv- 
port some heat capacitance, then as the piston moves through 
the cylinder the actual rate of heat transfer from the bulk 
gas to the voundary layer should instantaneously differ from 
the rate of heat transfer from the boundary layer to the 


eyliander wall. 


a Z a, Z a5 eee? femce Neat : (12) 


capacitance in boundary laye 


I,C. Phase Difference in Heat Flux 


The rate of neat transfer and therefore the loss expe= 
rienced from the system depends neavily upvnon the temperature 
difference between the working gas and the cylinder wall and 


upon the properties of the working gas. 
ae. iT as - Tyarl? (13) 
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During both the compression and the expansion stroke of a 
cycle, the system will pass through a point where the gas 
Pem@eravure woke equal the cylinder wall temperature. At 
Pmese two points the rate of heat transfer will equal zero 
as the temperature difference will equal zero ras eee 
= 0). Consider a system in which the cylinder wall boundary 
layer has zero heat capacitance. Here the temvoerature of 
the gas in the layer adjacent to the wall will be approxi- 
mately equal to tne bulk gas temperature. Therefore at tne 
points in the cycle where tne oulk gas temperature equals 
the cylinder wall temverature, the vorocess will ote adiabatic 
Emaecne race Of heat transfer will equal zero. 

However, if the system is assumed to have finite heat 
capacitance in the thermal boundary layer, the temverature 
Mumene fac in the layer adjacent to the wall could differ 
significantly from the bulk gas temperature. In this system 
[ifiemsOLnts Of zero rate of heat transfer will not occur 
where Se gas = Tall, ikee Vat eeeceur=ac the oie where 


Ty, aL 


wall. 
To illustrate this difference, if the compression ex-~ 


oundary layer 


pansion process is plotted on a pressure vs. volume diagram, 
Pie pOlNts Of Zero rate of heat transfer will occur where 
the orocess is adiabatic at that voint. Previously in this 
section it was shown that the process will be adiabatic at 
the voint where the trace is Coste Line where py > = 


constant. figure 4 illustrates a pressure vs. volume 


ne 





diagram for a system where tne points of tangency to the 
adiabatic lines coincide with the points where the bulk gas 
temperature equals the cylinder wall temperature. 

If a system with finite heat capacitance in the 
boundary layer is considered, the points at which the bulk 
gas temperature equal the cylinder wall temperature will 
not coincide with the points at whicn the system trace is 
tangent to the adiadatic lines. Figure 5 illustrates tnis 
difference. In both the compression and expansion stroke, 
a phase difference is identified that is the difference (in 
time or degrees of rotation) between the noints of actual 
zero rate of heat transfer in the system (voints where the 
process is verfectly adiabatic where  poundary layer = 
Teall? and the points of apparent zero rate of heat transfer 
in the system (voints where Ty 4. paces Twall)* 

Therefore the effect of the boundary layer on the heat 
monster in the cylinder must be determined in order to 
evaluate the imvact of the boundary layer on the total 
system loss. 

Other methods are certainly available that would facil- 
meace Calculation of the net cyclic loss and any phase dif- 
ference experienced in the points of zero rate of heat 
transfer. One of the most direct and accurate methods re- 
quires the plotting of a temverature vs. entropny diagram 


mom the cycle, Figure 6 is an example of such a trace. The 


area enclosed by the curve represents the net work done on 
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the working fluid which is also equivalent to the total loss 
Beapemtencea Guring tne cycle. For any single cycle the 


following statement should be true: 
prav =f Tds = loss. is 


On this temperature vs. entrovy diagram, the points on the 
trace that are tangent to a vertical line signify points 
where the change in entrovy is equal to zero. At these 
voints the vrocess is verfectly adiabatic and signify soints 
of zero rate of heat transfer. The temperature of the gas 
adjacent to the cylinder wall at these pvoints should be 
equal to the temperature of the cylinder walls. If these 
Permits Of actual zero rave of heat transfer differ from the 
points on the graph of apparent zero rate of heat transfer 


(where T = Ti .j,/> then some phase difference exists. 


bulk gas 
In Figure 6 a vhase difference is evident in both the com- 
pression and expansion strokes. 

The existence of this phase difference must be proven 
and then evaluated in order to provide information necessary 


to quantify the impact of the thermal boundary layer effects 


on the total system loss. 
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TT,4. xverimental Expander 


The apparatus used to investigate the theoretical pro- 
cesses of the expander was composed primarily of the motor 
and cylinder body of a modified air compressor. The active 
cylinder was bolted directly on top of the original com- 
meressor cylinder. This cylinder included tne gas charging 
port, temperature and pressure measuring instruments and nad 
maemadditional capability of water cooling the cylinder. A 
diagram of this apparatus is shown in Figure 7. The piston 
utilized to compress and expand tne gas was noused in the 
Meeiye Cylinder, The active cylinder bore was 2,00 inches 
and the piston stroke was 3.00 inches. During operation the 
piston remained within the cylinder, while its maximum up- 
ward travel at the top dead center vosition moved the piston 
Memtice Upper lin of the cylinder. 

B@ewcylianderm excension housed a stationary viston that 
formed the head of the cylinder volume. A threaded bolt 
extending through the top of the extension cylinder and into 
the stationary piston provided the additional capability of 
changing the volume ratio. By raising or lowering the sta- 
tionary piston through its four inch range with this 
threaded bolt, the cylinder provided volume ratios ranging 
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T1,3, Pressure ileasurement 


A strain gage oressure transcucer was utilized to 
measure cylinder pressure. The oressure measurement is 
determined in this type of transducer by the fact that the 
resistance of the transducer strain gages on the vressure 
sensitive surface changes in proportion to the woressure on 
the surface. The transducer output is then poronortional to 
the cylinder pressure. The transducer mounted in the cylin- 
der nead was of the range O to 100 posi and vas calibrated 
aes psi/mallivyoly for a 6.0 volt input. In order to 
Galculate the cylinder oressure, the transducer output at a 
fewmeoressure nad to be determined. Therefore vrior to 
meeconains any exnerimental data, the transducer input av 
atmosoheric oressure was recorded to provide the required 
oressure reference. To make this step easier a balance Dox 
was installed with the pressure transducer which vermitted 


Bae setting of a known pressure directly on the oscilloscove. 


T1.C. Volume Measurement 

The method utilized to measure the cylinder volume 
consisted of a circular cam mounted eccentrically on the 
compressor crankshaft. Because each cycle of the comvressor 
piston corresponded to one revolution of the comoressor 
crankshaft, the angle of a given spot on the circumference 
@eeune Shait was directly proportional to the position of 


the compressor voiston and therefore was also wvrovortional to 


a? 





the cylinder volume. A follower was mounted on the outer 
diameter of the circular cam so that the disvlacement of the 
follower was provortional to tne disvlacement of the viston. 
This system is displayed in Figure 8. 

hy er 7PCE displacement transducer was connected to the 
mollovwer and provided a D.C. ouUtDUT That las provortional to 
Biewdisplacement Of the follower, The calibration of the 
transducer provided the maximum output signal when the cylin-~ 
der volume was maximum at the piston vcottom dead center 
position. By knowing the maximum and minimum cylinder volume 
for the experimental run, the intermediate values for cylin- 
der volume could easily be calculated by using the transducer 


SUvouc. 


in. Cylinder ‘lall Temnerature Measurement 





Four thermocouvles were utilized to measure the average 
cylinder wall temperature during testing. Three thermocou-e 
ples monitored the temperature at different locations on the 
outside of the cylinder wall. The remaining thermocounle 
Mas placed on the outside of the cylinder cooling water 
jacket. The output of the thermocouples was displayed using 
an Autodata Nine Digital data logger that provided a contin-~ 


uous readout in degrees Centigrade. 


Mie. Dieital Storage Oscilloscove 


The outout of the pressure and volume transducers was 


provided as input to a Nicolet Digital Oscilloscove Model 
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206. The signal information vrovided to the scope was con- 
verted to digital form, stored in a buffer memory and then 
was disvlayed granhically. Digital information for each 
point remained available in storage. The Nicolet scope has 
the capability of storing up to 4096 data points. S5ecause 
both pressure and volume were being monitored, each storage 
cycle consisted of 2048 vressure points and 2048 volume 
Bomics., the scope had the capability of graphic disvlay of 
either the vressure or volume trace vs. time or an 4A-Y axis 
plot of pressure vs. volume. The time interval between data 
meiivus COULd be chosen on the scone to ensure that an entire 


cycle of data jould ve available in a single storage cell. 


Mier. Testing Procedure 
The following steps were verformed in this order during 
mee Gata acquisition run: 


1. The position of the stationary piston was set 
oroperly. 


2. The piston was moved to bottom dead center 
position and then checked by assuring that 
the output of the volume transducer was at a 
maximum 


3. The cylinder was ovened to the atmosphere and 
the vressure transducer outnut was set on the 
Ooscilloscove and recorded. 


4, The cylinder was charged with gas and the 
pressure transducer outout recorded. At this 
voint must ensure that the cylinder is not 
charged to a pressure that will over-vressurize 
the pressure transducer during compression. 


de the compressor motor was turned on and the 
oscilloscope set to automatically monitor a 
Single cycle of the cylinder. 
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The cylinder wall temperatures were recorded 
immediately upon comoletion of the sample 
cyclee 


The compressor motor was secured. 
The vressure vse volume data was monitored on 


the oscilloscope to ensure that the trace was 
satisfactorye , 
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Cie nett 
COMPULAR Re DUCBHON AMD INTERPRETATION OF DATA 





In order to efficiently interpret and display the raw 
data extracted from the Nicolet oscilloscope, the Institute's 
CMS computer system was used. The data was entered into oe 
computer files and then was operated on to provide various 
thermodynamic parameters at each point. Using this infor- 
mation thermodynamic traces for the cycle were plotted. 
Finally the losses for the cycle were calculated by numer- 


ical integration of the curves. 


moe A. Computer Invouts 
The digital pressure vs. volume data provided by the 
fmeobet oscilloscope was input into a data file that was 
accessible to the main vrogram. Other initial inputs re- 
quired included: 
Pines type Of veriinse gas used. 
Pee laewspeciile @as constant for the gas (R). 


3. The specific heat at constant pressure for 
the gas (C,)- 


4. The specific heat at constant volume for the 
gas (C.). 


De The cylinder volume at the tov dead center 
BOcmeron CVEDG) < 


6. The cylinder volume at the bottom dead center 
position (VBDC). 


7e The total number of data points (P, V pairs) 
entered into the computer data file. 


8. The time interval between data voints. 
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9. The olston revolutions per minute, 


ie tome cratunme readings from the cylinder 
wall thermocouvles in degrees Centisrade. 


Mmmewoe Calculavion of Reference Points 


In order to have the capability of accurately comparing 
data at varying volume ratios, compressor sveeds, and for 
different gases, computations performed on the data were done 
using common reference points. Therefore a standard method 
was required to determine each of the required temrerature, 
volume and voressure reference points. 

Cylinder oulk gas temperature was the first reference 
point determined and was set at 100 Se0 Ears temperature was 
chosen as it was a tynical average cyclic temperature experi- 
enced in ovrevious experimentation done with tne same exveri- 
mental cylinder apparatus. ‘nen exoressed on tne absolute 
scale the reference temperature was then set at the following 


value: 
T= Soom 8 


The method used to determine the reference volume con- 
Sisted of calculating the geometric mean using the volume of 
the cylinder at the top dead center and bottom dead center 
Peeemblons, The actual formula for the determination of 


reterence volume was 


a 
Ween 
Vy = (Re VTDC. (15) 
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Because the cylinder volume was expressed in units of cubic 
inches, the reference volume was also in cuodic inches. 

The determination of the reference pressure voint 
recuired the prior determination of the reference volume 
point. Using the reference volume point, the pressures 
corresponding to this volume for both the compression and 
expansion strokes were found. Tne reference pressure was 
then determined to be the average of the compression and 


expansion voressure values at the reference volume. 


> comp . Paxn 
a (16) 


on Za 
Reference oressure was exoressed in units of pounds force 


per square inch. 


ieee, oysceem riass Calculation 


In order to calculate the bulk gas temperature and 
entropy at each data point the system mass had to be deter- 
Mined. This calculation was performed using the ideal gas 
equation or state andthe reference temperature, volume, and 


pressure values. 
n= (17) 


The system mass was then expressed in units of vounds mass. 


meme. Calculation of Temnerature and sntrovv 


The bulk gas temperature at each data voint was then 


be 





determined utilizing the ideal gas equation of state and 


other previously determined parameters, 





C138) 


The temperature values were expressed on the absolute scale 
in degrees Rankine. 

The entrovy value for each data point was calculated 
referenced to the previously determined values of reference 
volume and pressure by the relation 


P. Va 
S. - Sy =~ AS= mC. ln| ==! + nC, val v, 2 (19) 


The entropy values were expressed in units of inch-pounds 


force per degree Rankine. 


Mite. Calculation of System Losses 

In order to calculate the system losses during the 
cycle, cyclic traces of pressure vs. volume and temperature 
VSe entropy were plotted. The area enclosed within either 


of these closed curves should equal the total system loss. 


Pay = Pas = loss (14) 


The integration was performed numerically using the trap- 
ezoid rule at each successive data point. Total area under 
the expansion curve was considered negative. The total 


system loss was then described in the units of negative 


b> 





inch-pounds force. 

An attempt was made to calculate tne area using a much 
more involved and typically more accurate five point Gaus- 
Sian quadrature integration method. However, vecause of 
the scatter of the digital data points near the ends of the 
curves, this method proved no more accurate than the much 


simpler trapezoid integration method. 


TiT,*. Non=-Dimensional Parameter Celculations 
PEm@icemucOncrOvilce the additional canability of com- 
Daring the computer output to the output data of any otner 
Petar apparatus, notwithstanding diiferences in the units 
involved, non-dimensional vcarameters were introduced. The 
computer routine tnen also included the calculations re- 
quired to vroduce the outout of non-dimensional values of 
oressure, volume, temperature and entrovy. Plots were also 
made of non-dimensional oressure vs. volume and non-dimen- 


Slonal temperature vs. entropy in order to determine the 


non-dimensional loss values. 


TitT,Y.1, Non-Dimensional Volume 
The volume data voints were nonedimensionalized by 
dividing each data point by the total displaced volume of 


the cylinder. 


v 
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iietec. JOl—Dimensional Pressure 

The cylinder pressure data points were non-dimension- 
alized by dividing eacn data noint by a mean oressure cal- 
culated from the comoression stroke. This mean pressure was 


defined as the total work added during tne compression 


stroke divided by the total disolaced volume of the cylinder. 


he 
= “comp 
Pep = (HBDCAVTDC) (21) 


wach oressure data point was then divided by this mean ores- 


sure during comoression to find the non-dimensional oressure 
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a oe = (22) 


Mme. so. NOn—Dimensional Temoerature 





The bulk gas temperature data voints were non-dimension- 
alized by dividing each data point by the reference temper- 
ature. 


. 
TSTR, = C25) 


a 
O 
Tii.*.4. Non-Dimensional =ntrony 
The entropy values were non-dimensionalized by dividing 


the calculated entropvy values by the mean pressure during 


compression and the displaced volume of the cylinder and 
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PecnemuLctiolyine by the reference Cemperature. 


S5 Th 
SSTR = DP UUBDCLUTDC) (el) 
Pen /BDC=-VTDC 
Non=-Dimensional System Losses 





Nonedimensional values of total system loss were cal- 
culated in the same manner as the dimensional loss values 
were determined. Flots were made of non-dimensional vores- 
sure vse volume and non-dimensional temperatures vs. entrony 
and the enclosed areas were calculated by numerical inte- 


em@avtOn utilizing the trapezoid method. 


mec. COMouter Outour 

A listing of tne Fortran comouter program is displayed 
maepoendix A. A listing of the test data file used in the 
computations is presented in Appvendix 8B. 

The output of the computer routine includes the fol- 
lowing data: 


1.e Pressure, volume, bulk gas temperature and 
SMmerewy ac each data cvoOLnt. 


ee Non-dimensional pressure, volume, bulk gas 
temperature and entrovy at each data point. 


5e Piston revolutions ver minute. 
4. Time interval between data points. 


5e Reference temperature volume and pressure 
values. 


6. Average cylinder wall temperature. 


7e System mass. 
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Pressure vse volume vlot. 

a : , 
Temperature vse volume plot. 
Temperature vs. entrovy nlot. 


Non-dimensional pressure vs. non-dimensional 
volume plot. 


Nonedimensional temoverature vs. non-dimen- 
Sional volume ovolot. 


iion<dimensional temperature vs. non-dimen- 
sional entropy plot. 


Seleculacaems Of camlensional and non=dimen- 
sional system logs. 


in example of the comouter output is disolayed in Anvendices 


C 


and D. 
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The results of the single data set used are oresented 
in the computer outout in Appendices € and D. Helium was 
utilized in the cylinder as the working gas for the exver- 
iment. 

The final oressure vs. volume and non-dimensional 
pressure vs, non-dimensional volume plots resulted in 
extremely smootn traces. tach depicted what could easily 
be identified as a conventional pressure vs. volume trace 
for this tyoe of apparatus. Basically the same comment 
could be made about both the temperature vs. volume and 
non-dimensional temperature vs. non-dimensional volume 
vlots. However, a greater degree of scatter seemed to be 
evident in the temperature vs. volume vlots than in the 
oressure vse volume plots. 

In the temperature vs. entroovy and non-dimensional 
temverature vs. non-dimensional entrovy plots tne scatter 
is much more evident. The shave of the trace also is 
somewhat different from a shape that one might expect the 
plot to have after first seeing the vressure vs. volume and 
temperature vs. volume traces, Figure 9 is a single line 
Swve aporoximation of the TSTR vs. SSTR computer trace 
from Appendix D. The figure shows the adiabatic and iso- 
Piermal limits of a plot of this tyne. The adiabatic limit 


Deing a curve where the comoression and exvansion strokes 
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approach being perfectly adiabatic while the isotnermal 
limit aporoaches a curve where the compression and expnansion 
Surokes are veriectly isothermal. The [StR vs. SSTR volot 
has been dratm on a gravh with axes of eaual lengtn so that 
tne curve nas not been mistakenly lengthened in either 
Girection to give the false appearance of deing more acila- 
batic or isothermal. Nevertheless, the trace seems to 


approach the adiabatic limit in the high temperature portion 


iy 


@eecne plot, «nile approaching the isothermal limit in tne 
low temperature vortion of the ovlot. On the comcression 
stroke tne trace becomes verfectly adiaoatic at a point 

7 : ; TMCanD «= oe a =u 1 
~oimm eaa) tne seroks at TSTR = .665 (T = 48h-R). On the 


expansion stroke the trace also becomes perfectly adiabatic 


O Oo 
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mea OOint carly in the stroke at TSTR = 1.15 (T = 644 
At each of these two voints that actual rate of heat 
transfer in the system through the cylinder walls equals 
ZCLO. 

The points of apparent zero rate of heat transfer occur 
at T= aa = Eeone Re The difference between these points 
and the points of actual zero rate of heat transfer is the 
ohase difference. In both the compression and expansion 
strokes the ohase difference is leading as the noint of 
Beever Zero rate of heat transfer occurs nrior to reaching 


the voint of anparent zero rate of heat transfer. The 


actual phase differences can be comvuted to be: 


LO 





- , O ‘ 
compression stroke a) leading 


expansion stroke 64° leading 
The actual system loss figures comouted were: 


PdV loss = #24.931 in-lo. 


TdS loss = =25.107 in-lo. 
Non=dimensional PdaV loss = -0.0584 


Non-edimensional TdS loss = ~0.0588 


in both the dimensional and non-dimensional cases, tne PdV 


and TdS losses showed less tnan one vercent difference. 
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fee Data Point Scatter 

The scatter evident in the temperature vs. volume plots 
and particularly in the temperature vs. entrovy plots 
deserves explanation. In these plots this scatter occurs 
primarily in the center section of the trace, away from the 
areas Of minimum and maximum cylinder volume. At the unper 
and lower end vortions of all the plots, the data soints are 
Senerally packed tightly together in an orderly fasnion, 
Mies is due co the fact that at the ends of the compression 
and expansion stroke, the piston is decelerating rapidly, 
nalting abruptly and then again rapidly accelerating away 
mom 1tbs extreme position. This stopping and starting 
motion requires a great deal of time to cover a small cis- 
tance in comparison to the viston's much more ravid vace 
at the mid-stroke position of its travel. Because tne data 
points are taken at equal time intervals, many more data 
points are taken at the extreme positions of the plots 
providing the tighter packing and the generally smoother 
apovearance in those areas. 

Near the mid-stroke positions of the cycle a different 
situation exists. Here the piston is moving much more 
Maplaly relative to the velocity experienced at the extreme 
cylinder vositions. Again as the data points are taken at 


equal time intervals, a greater distance should exist 
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between points on each trace near the mid-stroke position. 
Pies Situetion alone does mot account for the scatter, but 
along with two other possible causes, the existence of tne 
scatter can be better explained. 

First, the Nicolet oscilloscope samovles 2043 volume 
and 2948 pressure data points per scope storage cycle at 
equal time intervals between each point. The procedure of 
the scove is to sample one volume data point, wait the 
aporopriate time interval and then samole one vnressure data 
point, continuing in this manner until the storage cycle is 
full. The scove does not nave the cavability to sample 
pressure and volume simultaneously. Therefore tne resulting 
oressureevolume data point pair formed by matching a volume 
Gata point with the following oressure data point does not 
Beeculy cCOincide with that voressure-volume data point vair 
that would occur if the cavability existed to sample each 
pair simultaneously. This fact brings a small amount of 
Bemor INTO the position of each data voint. 

Second, the digital nature of tne data generates addi- 
mammal error. The accuracy of the digital data outout of 
Pm@e Nicolet was limited to .05 millivolt. ‘Jith the large 
number of data points taken, the changes in outout voltage 
between many points was less than .05 millivolt. Therefore 
voltage changes of less than .05 millivolt were rounded off 
to the next .O5 millivolt increment creating some error in 


the data. This might have not been critical except for the 
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fact that the volume data was being read on the 100 milli- 
volt scale on tne oscilloscope, and the range of the data 
covered only 3.00 - 17.90 millivolts (roughly 15% of the 
range of the scale). The 100 millivolt scale was the best 
scale to use as it was the smallest scale available on the 
scope that would accept the range of the volume data. Howe 
ever, a 4X or 5X signal amplifier should have been connected 
mm the volume signal wiring to utilize the full range of the 
100 millivolt scale. This would nave alleviated the problem 
epeesetns due tO the .O5 millivolt accuracy of the Nicolet 
oscilloscotve. 

Migsem lini tavrtons Of the disital oscilloscove in con- 
mamevlon wath the magnitude of the midestroke piston velocity 
producing a greater distance between data voints in this 
area Of the vnlots, are surely a sufficient basis with which 


wemccetemor tO justify the scatter, 


V.8. Shane of Temverature vs. Entrovy Plots 


The actual shape of the temperature vs. entropy plots 
exhibits a few variations from what might be expected after 
seeing the other plots. In the high temperature end of the 
plot, the curve seems to approach its adiabatic limit, 
especially on the initial portion of the expansion stroke 
(Figure 10). The explanation for this is fairly simple. 

At the start of the exvansion stroke, the temverature of the 
bulk gas and the temperature of the gas in the boundary 


layer are both greater than the cylinder wall temperature. 
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However, the temverature of the bulk gas is still signifi- 
eanvingeomeaver Gian the boundary layer gas because the lover 
memocracure cylinder wall has a greater effect on the temper- 
ature of the gas in the thin boundary layer. As the expan- 
sion stroke begins, the vressure drops and the bulk gas 
temperature drops rapidly. The temperature of the boundary 
layer gas drops also, but it drovs at a much slower rate as 
it started the expansion stroke at a significantly lower 
temperature. Because the temperature difference between the 


boundary layer gas and the cylinder wall is not great for 


Kt 


mie tirst thirty degrees of tne expansion stroke, very Little 
neat is transferred out of tne system and the pvnrocess av-= 
proaches its adiabatic limit. 

At tne low temperature end of tne plot, near the end 
of the expansion stroke the curve seems to apvoroach its 
moetnermal limit (rigure 10). The explanation fOr cds aS 
yery Similar to that for the near adiabatic case. Near the 
end of the expansion stroke the bdoundary laver gas temper- 
ature has fallen somewhat lower than the cylinder wall tem- 
perature. In the last thirty degrees of the e&pvansion stroke 
the volume in the cylinder increases, but it increases at a 
relative rate much less than the rate of increase near the 
start of the expansion stroke. Therefore the change in the 
oulk gas temperature will also change at a rate much less 


than the rate of change at the start of the expansion stroke. 


Even though the process aporoaches its isothermal limit, 


46 





there is still a sufficient temperature difference cetween 
the boundary layer gas and cylinder wall to pvermit signifi- 
eant heat transfer from tne system. 

The temperature vse entrooy plot also has what seems to 
be an inflection near the rete eele position on the compres-~ 
sion stroke (Figure 10). At this inflection there is an 
abrupt change in the slone of the compression trace. Because 
of its proximity to the mid-stroke position on the curve, 
mms inflection may be due to the change of viston accelera- 
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tion to deceleration at midestroke. Wo such ¢g 
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men 1s evident in tne exmansion curve. 


weer loam iteeet ion gf coe £clieoatic Points in the Crecle 





miecmraaentiiicatzon Of the tvo points in the cycle that 


are adiabatic is extremely important in determining the 


cr 


1 
ie 


actual ohnase difference vetween the heat transfer and 
bulk gas to wall temperature difference. It is exvected that 


both the compression and expansion strokes should each have 


es 


: : : vy T 
Saeeadiabatic point where boundary layer wall, 


Ideally 
this should occur at the voints of maximum and minimum values 
of entrovy. However, odecause of the scatter in the tempera- 
fee TS. entropy nlots it was difficult to envision the 
actual location where the curve became tangent to the ver- 
jeetie inerefore to identify the location a rough curve was 
passed through the data points (Figure 9) and the adiabatic 


points were then evaluated at the locations where the curve 


Massed through the vertical, 


4’? 





; , nat eS 
Giemerccuswenmesewaoctaneatic locations, Lines or PV = 


constant (adiabatics) were constructed on the pressure vs. 
volume plot (Figure 11). The locations on the trace that 
were tangent to an adiabatic line vould signify an adiabatic 
feet. Lae results of this procedure verified the selection 


™ 
4 
— 


eeceae adiabatic voints in Figure 9. 


V.D. Effects of the Boundary Laver 


The effects of the boundary layer on the rate of heat 
transfer, or more specifically on the losses in the system 
Beemmost important. It is obvious that it is the temperature 
of tne boundary layer gas and not tne dulk gas temperature 
meen 1s tae driving force on the transfer of neat through 
Giemcylinder wall, ii the boundary layer gas temperature is 
Sreater than the temperature of the cylinder wall, the system 
Will lose heat. If the boundary layer gas temperature is 
less than the temperature of the cylinder wall, the system 
will gain neat. 

Along the compression stroke and part of the expansion 
stroke from the minimum to maximum entropy value the boundary 
layer acts to cool the system (Figure 10). On this area of 
the curve the boundary layer gas temperature is basically 
Sreater than the cylinder wall temperature and heat will be 
transferred out of the system through the cylinder walls. 
M@emacat transier in this portion of the curve is driven by 
the change in volume and not the change in the bulk gas ten- 


perature of tne system. 
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Along the expansion stroke from the minimum entroovy 
value to the maximum entrovy value the boundary layer acts 
to heat the system (Figure 10). The temperature of the bulk 
eas and the boundary layer is falling; however, the voundary 
layer temverature is generally less than the cylinder wall . 
temperature and therefore heat will be transferred into the 


system through the cylinder walls. 


Vea. sxolanation of the Phase Difference 

The ohase difference that exists in the zero rate of 
neat transfer location is a very important effect of the 
poundary layer (Figure 9). The explanation of the cause for 
Mee DJnase diiference should ove fully discussed. 

Tne oulk gas temperature is an imvrortant sarameter in 
this tyne of system. Secause the bulk gas temvnerature is 
much easier to measure and to understand than the more dif- 
ficult to measure boundary layer temvoerature, the oulk gas 
memperature iS utilized much more in this tyve of system, 
Therefore, if the bulk gas temperature is utilized, then 
when deriving the apvarent expression that troould define the 
rate of heat transfer from the system, the temperature ex-= 


tremes utilized would be the bulk gas temverature and the 


cylinder wall temperature and the expression would be: 


Q e 


ee (Toulk Bas Mwall?* 


memeonisS expression the points of zero rate of heat transfer 


are the coints on the cycle where T, 
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The existence of a boundary layer with finite neat 
S-wacttvance alvers the former incverpretation of the problen, 
Mae temperature of the gas in the layer next to the cylinder 
Wall will not equal the bulk gas temperature, out will eaual 
the boundary layer gas temperature. Tnerefore, the expres- 
sion defining the rate of heat transfer to the cylinder walls 


Wa becone: 


= x 7 oo a a 
a = \T poundary layer one 


This exvression svecifies that tne noints of zero rate of 


neat transfer would oe the noints on the cycle where 


m 


mM hy 1 + ’ i 2 ” 
"poundary Layer 7 lwail, Because the oulxK gas temperature 


is not eaual to the boundary layer gas temnerature in this 
system, some caifference till exist in the time of apparent 
zero rate of heat transfer and actual zero rate of heat 
cransfer. This difference is tne zero rate of heat transfer 
ohase difference. ‘The pnase difference in both the com- 
oression and expansion strokes of this experimentation was 
leading, or the actual zero rate of neat transfer occurred 
in the cycle before the apparent zero rate of heat transfer 
(Figure 9). 

The explanation for the phase difference exists in the 
examination of the temperature history of the boundary layer, 
As the compression stroke is started the boundary layer gas 
temperature is slightly greater than the bulk gas temperature. 


Mis is due to the fact that the cylinder wall temperature is 


5 





greater than the bulk gas temperature and the cylinder wall 
temperature nas a much greater effect on the temnerature of 
the boundary layer gas rather than the bulk gas. As the gas 
is comcvressed the bulk gas and boundary layer gas temperatures 
rise, but because the vcoundary layer gas temperature started 
the stroke at a hisher temperature than the oulk gas, the 
boundary layer gas temperature passes through the cylinder 
wall temperature vrior to the bulk gas oroviding a leading 
onase difference. 

As the expansion stroke is started botn the bulk sas 
Pye oOoundary layer gas are cooled oy the drov in voressure. 
Sut ovecause the ooundary layer gas temverature vas signifi- 
cantly lower than the obulk gas temperature at the beginning 
of the exvansion stroke, the boundary layer gas temperature 
vasses through the cylinder wall temperature vorior to the 


fe 


bulk gas oroviding a leading ohase difference. 
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CON CeUSrOus 


One of the original objectives of the work included 
the capability of direct transmission of the data from the 
experimental exvander to the Institute's CMS computer system. 
This canability required tne use of a converter that could 
vroverly interface with the Nicolet storage oscilloscove and 
the CMS comouter system, The zlectric Povrer Systems Labcora- 
tory was assigned the task of designing and assemoling this 
mmeeriace unite Due to mannower shortages, tne zlectric 
Power Systems Laboratory was not able to comclete the as- 


memory Of ths interlace wunit orior to the suomission of tnais 


by 


vaper. the unavailability of the interface unit required the 
Pewee tO De read point oy point digitally from tne Nicolet 
storage oscilloscove and then manually entered into a CHS 
computer system data file. 

Behe bae Use Of an Operable interface unit, dozens of 
cycle runs could have been performed. Eowever, the lack of 
the interface significantly limited the number of cycle runs 
that could be evaluated. Further setoacks occurred when 
serious oroblems were encountered in the digital and graphic 
data presentation from tne Nicolet storage oscilloscone. The 
pedal resuit was the availability of only one good set of 
cycle data for analysis. 

Although the technique develoved in this vaper to 


analyze the data was used on only one set of experimental 


Do 





data, the results obtained oroved that the analysis method 


ry 


was highly credible. ‘Inetner one or one hundred sets of 
cycle data nad ocoeen evaluated, the same method would nave 
been utilized. 

The additional runs that could have been verformed ‘vith 
an overable interface unit would have greatly increased the 
scope of the experiment. Parameters sucn as piston RPM, 
Reynold's number, cylinder volume ratio, cylinder vressure 
and type of working gas could have oeen varied. The results 
that would have ceen voresented with the variations in the 
above varameters would have orovided insight into the actual 
meet transtem losses and ohase difference for a wide range 
Me conmcivions, This general information could then be 
extravolated for use in the design of other devices of this 
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